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ABSTRACT 

Data products from the Advanced Camera for Surveys Virgo Cluster Survey 
are used to understand the bulge star formation history in early-type galaxies at 
redshifts z > 2. A new technique is developed whereby observed high-redshift age- 
metallicity relationships are utilized to constrain the typical formation epochs 
of metal-rich or "bulge" globular clusters. This analysis supports a model where 
massive Virgo galaxies underwent an extremely intense mode of bulge globu- 
lar cluster formation at z ~ 3.5 that was followed by an era of significant bulge 
growth and little globular cluster production. Intermediate-mass galaxies showed 
a less-intense period of globular cluster formation at z ~ 2.5 that was synchro- 
nized with the bulk of bulge star growth. The transition between the massive 
and intermediate-mass galaxy star formation modes occurs at a galaxy stellar 
mass of M s teiiar ~ 3 x 10 10 Mq, the mass where many other galaxy properties 
are observed to change. Dwarf early-type galaxies in Virgo may have experienced 
no significant period of bulge globular cluster formation, thus the intense star 
bursts associated with globular cluster formation may be difficult to directly ob- 
serve at redshifts z < 4. Though the above conclusions are preliminary because 
they are based upon uncertain relationships between age and metallicity, the 
technique employed will yield more stringent constraints as high-redshift galaxy 
observations and theoretical models improve. 

Key words: galaxies: bulges < Galaxies, galaxies: clusters: individual:... < 
Galaxies, galaxies: formation < Galaxies, galaxies: high-redshift < Galaxies, 
galaxies: star clusters < Galaxies 



1 INTRODUCTION 

Globular star cluster formati on only occurs during 
strong star formation even t s (lHarris fc Pudrita Il994l ; 
lElmegreen fc Efremovl 1 19971 : iLarsenl 120091 ). Most GCs 
found in earl y- type galaxies formed a t redshifts z > 2 (see 
references in iBrodie fc Str adcr 2006) , thus globular clus- 
ter (GC) systems are valuable observational tools to help 
understand the nature of major star formation events in 
early Universe dAshman fc Zepflll992l:lForbes et al.lll997l ; 



ICote et al.lll998l : lHarrisll200ll : IBrodie fc Straderll2006h . 

High-redshift galaxy observations continue to pro- 
vide more and more detailed galaxy properties (e.g 
ages, metallicities, sizes, star-f ormation rates) at the red- 
shifts of GCs f o rmati o n (e.g. Hopkins fc Beacoml 120061 
iBouwens et~ai1 120071 ; iReddv et all 120081 : iFranx et al 
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20031: Ivan Dokkum et ail 120031; iGlazebrook et al.l |2004|: 



Daddi et al.l 120051: Icimatti et all 120081 : llvlaiolino et al.1 
20081 : van Dokkum et alT 2010t ). Though much progress 



has been made in this field, interpretations are gener- 
ally limited due to the challenging nature of such ob- 
servations. By combining constraints from near-field GC 
system observations with far-field galaxy observations, 
certain limitations can be overcome and unique predic- 
tions for this i mportant period of galaxy formation can 
be made (e.g. IShapiro et alj|20ich . 

One of the strongest predictions resulting from GC 
system work is that many galaxies experienced two modes 
or epochs of intense star formation sometime before z ~ 
2. This follows from observations that early-type galaxies 
with stellar masses of M s t e iiar > 10 10 Mq generally host 
two GC metallicity su bpopulations (e.g. IStrader et al.l 
120061 : iPeng et al.ll2006l ). The metal-poor or "halo" GC 
subpopulation is thought to have formed within the 
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numerous metal-poor proto-galaxies that likely domi- 
nated the early Universe. The metal-rich or "bulge" 
GC subpopulation likely formed later and shows similar 
properties to the ho st galaxy's bulge (e.g. metallicities , 
spatial distributions [Harrid l200ll: IForbes fc Fortell200ll; 



Harris fc Harris! |2002|; IDirsch et al.l 120051 ; IBassino et all 
20061 : iForte et al.ll2007l . 120091 ). 

Thus by directly comparing the properties of a GC 
subpopulation to the galactic component that they are 
associated with, the relative contribution of e.g. bulge 
stellar mass produced during the intense MR GC forma- 
tion events can be understood. Furthermore, this type of 
comparison can constrain more general galaxy formation 
models by also considering observational and theoretical 
analysis of high-redshift galaxies. 

The present work conducts such an analysis and 
compares properties of bulge GC subpopulations to their 
host galaxy's bulge stars. A model of bulge star for- 
mation is developed for redshifts z > 2, which is com- 
pared to existing constraints from direct, high-redshift 
galaxy observations. The GC system observations come 
from a imaging survey of 100 early-type galaxies in the 
nearby (~ 15Mpc) Virgo Galaxy Cluster using the Ad- 
vance Camera for Surveys mounted on the Hubble Space 
Telescope (HST). The sample size, its homogeneity and 
broad range of high-level data product resulting from 
the Advanc e Camera for Su rveys Virgo Cluster Survey 
(ACSVCS; ICote et all [20041 ) allow for a detailed inves- 
tigation into the relationship between galaxies and their 
GC systems, hence constraints on the first few gigayears 
of bulge star formation in Virgo early-type galaxies. 



1.1 Globular Cluster System Observations 

Spectroscopy work to age-date individual GCs in ellipti- 
cal galaxies has shown that the bulk of GCs are very 
old, with typical ages of > 10 Gyrs or z > 2 (e.g. 



Larsen et al.ll2002l: 


Strader et al. 2005: Puzia et all 


2005; 


Beaslev et al.l 2006 


:IConselicell2006 


; Sharina et al.l 


2006; 


Pierce et al. 2006 ?; ICenarro et al. 


2007; Beaslev et al. 


20081). Because GCs are dominated bv old stellar dod- 



ulations, the colour distribution of a galaxy's GC sys- 
tem can be in terpreted as its intrinsic metallicity dis- 
tribution (see IStrader et alJl2007l : iKundu fc Zepli2007l : 
ISpitler. Forbes fc Beaslev 20081 ). Observations show that 
most extrag alactic GC systems have a bimodal colour 
distribution (jStrader et al.ll2006l : IPeng et al.ll2006h . which 
implies they host two metallicity subpopulations: metal- 
poor (MP) GC subpopulations with typical metallicities 
of [m/H] r 2 to -1 and metal-rich (MR) GC subpop- 
ulations ranging [m/H] ~ — 1 to 0. 

MR GCs show roughly similar chemical proper- 
ties to the their host galaxy's bulge M ste iiar_ (e.g . 



Harrisll200ll; IForbes fc Fortell200ll ; lHarris fc Harrisll2002i : 
Forte et al.ll2007l . 2009) and the spatial distributions of 



MR G Cs resemble the hos t galaxy's stellar distribution 
(e.g. IBassino et al.l 120061 ). though the GCs generally 
show a shallower distribution in the central regions of the 
galax y possibly due to GC destruction (e.g. IDirsch et alj 
120051 ). These observations suggest the formation of MR 
GCs and galaxy bulges are closely linked. 
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Figure 1. (g — z)q co lour profiles of Virgo e arly-type galaxies 
from the Sersic fits of Fcrrarcsc et al. I teood) . For each galaxy, 
its mean MR GC colour is subtracted from the galaxy's colour 
profile to illustrate the relative enrichment of these two com- 
ponents. Top, middle and bottom pane ls correspond to galax- 
ics in the ACS Virgo Cluster Survey dCote et al.ll2004l ) with 
^stellar ^ 10 10 Mq (hereafter "low-mass" galaxies), 10 10 < 
M st eii a r < 10 n M Q ("intermediate") and M stMar > W n M Q 
( "massive" ) , respectively. Colour profiles with magenta dashed 
lines have central blue (young) cores and are removed from 
further analysis (see Ej2j . The ±1<t colour boundaries for a 
typical MR GC subpopulation in that galaxy mass bin are 
represented by the two dotted lines. All massive galaxies show 
colours that are redder (i.e. more enriched in metals) than the 
typical MR GC. Intermediate and low-mass galaxies show a 
range of differences. 



MP GCs s how metallicities m ore resembling galaxy 
halo stars (e.g. ISearle fc Zi nn 1978) and are thus thought 
to have formed during a "pre-galactic" era before galax- 
ies started to assembly in a ACDM universe. Because 
the intrinsic metallicity distribution of early-type galaxy 
stellar metal-poor halos have only been studied in a 
few galaxies (e.g. lElsonl 1997; lHarris fc Harris! 1200: 



Harris et all 12007 



2009; 



Foster et al 



lElsi 
iNor 



Norri s et al.l 
12Q09J), 



2008; Weiimans et al 



the relationship between MP 
GCs and their host galaxy is more difficult to discern 
and hence constraints from MP GC subpopulations will 
not be considered here. 



2 ACSVCS DATA SAMPLE 

The properties of MR GCs, associated with Virgo Galaxy 
Cluster early-type galaxies, are investigated in the follow- 
ing sections. A Virgo GC system is considered to host a 
MR GC subpopulation if the GC system's colour (hence 
metallicity) distribution yields can be we ll-represented 
by two Gaussians (see iPeng et al.l 120061 ) and the fi- 
nal MR GC number estimate must not be equivalent to 
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zero within the measurement and co ntamination errors . 
Thi s informatio n is co ntained within IPeng et al.l (|2006l ) 
and IPeng et al.l (|2008l ). To prevent contamination from 
nearby GC systems, galaxies are not considered if are em- 
bedd ed into the proje cted GC system of a larger galaxy 
(see IPeng et alj|2008h . GCs are mostly old (see © and 
contain little dust (e.g. iBarmbv et al.ll2009h . thus their 
colours are largely determined by their intrinsic metallic- 
ities. The ACSVCS GC colours (F475W-F850LP; here- 
after g — z) are converted into a metalliciti es estimate 
using the empirical relation given as eqn. 2 in lPeng et al.l 
(|2006l ). This is valid for an old stellar population and is 
based upon direct observations of Milky Way and M87 
GCs. Mean GC g — z col ours and measure d intrinsic 
colour dispersions are from lPeng et al.l (|2006h - 

Mean galaxy metallicities of the ACSVCS sample 
are critical for much of the following analysis. Com- 
pared to the simple stellar populations of GCs, ob- 
served galaxy light comes from multiple stellar compo- 
nent components and is thus more difficult to inter- 
pret. Identifying the population of stars that formed 
during the same epoch of GC formation and then un- 
derstanding their intrinsic metallicity content is ham- 
pered by the w ell-known age-metallicity degeneracy (e.g. 
IWorthevlll99i ). Spectroscopic analysis can help break 
this degeneracy, but such work is generally limited to 
the central regions of a galaxy and are currently avail- 
able for only a small subset of the ACSVCS sample. 
Spectroscopic observations have show n that early-type 
galaxies are generally as old as GCs dKuntschner et alj 



20021: iThomas et al.ll2005l: ISanchez-Blazauez et al.ll2006l : 
Koleva et all 120091 ; ISmith et aDl2009h . though interme- 



younger central stellar populations (e.g. Geha et al. 


2003 


;ISanchez-Blazc 


uez et al.ll2006l;ISansom & Northeast 


2008 


; Smith et al.l 


2009). However, such observations 



can be influenced by a relatively small (in terms of 
M s t e iiar ) young star burst that can briefly domi- 
nate a galaxy's central light (e.g . Lisker et al.l [2006; 
iTrager et all I2OO8I ; IProctor et aTT 120081 ). In" any case, 
a subset of the low-mass elliptical population shows 
ages that are identical to old massive galaxies (e.g. 



Kuntschner et al.ll2002l; IKoleva et alj|2009l ; [Mendel et all 
20091 ; iTolstov et alj|2009l ). 



Analysing galaxy photometry is one way to under- 
stand the stellar populations of the ACSVCS galaxies in a 
homogeneous fashion. Though galaxy broadband colours 
like g — z are prone to uncertain interpretation, they 
do probe the global properties of a galaxy and there- 
fore should be se nsitive to spatial variations in the stella r 
populations (e.g. IKoleva et al.l 2009; Spo laor et al.|[2009l ). 
which are sometimes missed in spatially-unresolved spec- 
troscopic analysis. In order to study the ACSVCS sample 
in a homogeneous fashion, a galaxy's (g — 2)0 colour is 
transformed into a metallicity estimate by assuming the 
stellar population is dominated by stars of a similar age 
to GCs and using the same colour-metallicity transfor- 
mation that is used for the GCs (see above). This trans- 
formation is inappropriate for galaxies that show young 
stellar populations in their photometry. Those with ob- 
vious young components (see below) are culled from the 
present analysis. Despite this precautionary action, it is 



possible that a subset of the remaining ACSVCS sam- 
ple could have luminosity- weighted ages that are younger 
than a typical GC. The metallicity values derived for such 
galaxies would therefore be lower limits. The reader will 
be r eminded of this possib ility when appropriate. 

iFerrarese et alj (|2006h identify two galaxies in their 
sample (VCC 798/M85 and VCC 1499) that show sig- 
nificant young stellar components. They also identify 
certain galaxies with corrupt or tidally truncated light 
profiles, which are also excluded here. A few of the re- 
maining galaxies in the IFerrarese et al.l (|2006l ) work have 
blue cores in their colour-profile, consistent with a br ight, 
young stellar population (e.g. iLisker et al.l l200fj ), as 
shown in Fig. [T] Since a bright, young core complicates 
the estimation of a galaxy's M s t e iiar and average metal- 
licity, these galaxies are excluded from the analysis. In 
particular, if a galaxy's global (g — z)o colour is 0.2 mag- 
nitudes smaller than their colour at 10 effective radii (or 
the outermost radial bin covered by the HST ACS field) 
they are excluded. Galaxy (g — z)p colour prof il es an d 
z— band effective radii are from IFerrarese et al.l (|2006l ). 
Fig. Q] shows the ACSVCS galaxy colour gradients sub- 
tracted by each g alaxy's mean MR GC colour taken from 
IPeng et alj i|200rj) . Those with blue cores are highlighted 
in the Figure. 

Galaxy stellar masses are from 2MASS K-band lu- 
minosities (which is largely insensitive to metallicity ef- 
fects), and a mass-to- light ratio (M/L x = 0.86) appropri- 
ate for an old stellar p opulation (see ISpitler et al. I l2008l ; 
ISpitler fc Forbesll2009l ). 



3 CHEMICAL DIFFERENCES BETWEEN 
MR GCS AND BULGE STARS 

As reviewed in Section ri.il MR GC subpopulations share 
some properties with their host galaxy's bulge stars. In 
this section, the chemical similarities of the MR GCs and 
galaxy bulges are investigated. 

3.1 Analysis 

In Fig. [T] the difference between the galaxy colour gra- 
dient and its mean MR GC subpopulation colour is pre- 
sented for ACSVCS galaxies. Massive galaxies show red- 
der colours than the typical MR GCs by ~ 0.1 mag, while 
lower mass galaxies show a spread of differences. In the 
same Figure, the average ±lcr colour dispersion of the 
MR GCs in each M ste uar bin are presented as dotted 
lines to illustrate the range of colours characteristic of 
the MR subpopulations (<r 9 _ z ± 0.15, ±0.14 and ±0.13 
mag. in order of the massive to low-mass galaxy bins). 
Galaxy colours generally overlap with the observed in- 
trinsic colour dispersion of their MR GC subpopulations, 
which is why previous work (see Hl.ip has claimed that 
bulge stars and MR GCs are roughly similar in terms of 
colour and hence metallicity. 

A more quantitative investigation of this relationship 
is presented in Fig. [2] where the average metallicities of 
both components are given. Galaxy metallicity reflects 
the luminosity-weighted global colour computed from in- 
tegrating the best-fit Sersic profiles in each band to infin- 
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Figure 2. Mean MR GC (open circles), MP GC (open squares) and galaxy (solid circles) metallicities against host galaxy 
stellar mass. At approximately M ste u ar ~ 10 10 Mq, the MR GC and galaxy mass-metallicity relationships start to merge. For 
M s t e iiar ^ 3 X 10 10 M 3te n ar , the two components are typically offset by a factor of ~ 2 in metallicity. Error bars are la- 
measurement uncertainties. Magenta points are the Milky Way data. 



ity and the assumption that the colours are dominated 
by an old stellar population (see JZ}. It is apparent in 
Fig. [5] that both components become more enriched with 
host M s teiiar , reflecting the well-known mass-metallicity 
relationship in gala xies (e.g. iKodama fc Arimotolll997l ; 
Gallazzi ct al. 2006) and an established trend among MR 



GCs (iForbes et all 19971 ; lLarsen et ai1l200ll ; IStrader et all 



l2006l : |Peng et al.ll2006h . 

Also included in the Figure are MP GC mean metal- 
licities. The difference between the sub population metal - 
licities remains relatively constant (e.g. iPeng et al.ll2006l ) 
until lower galaxy masses where MR GC subpopulations 
are no longer detected. It is still possible small numbers of 
MR GCs are present in these galaxies, but spectroscopic 
confirmation of the few MR GC candidates is required. 

From Fig. [2] it is apparent that the mean metallici- 
ties of the bulge stars and MR GCs are fairly well sepa- 
rated for massive galaxies, but become more similar for 



lower mass galaxies. For massive galaxies, the bulge stars 
show more enriched stellar populations than their typ- 
ical MR GCs by ~ 60% in metallicity. For galaxies of 
M stMav < 3 x 10 10 M Q , the average difference between 
the two components apparently decreases, as the MR GC 
and galaxy mass-metallicity relationships begin to con- 
verge. These galaxies have ~ 15 MR GCs on average, 
so the mean MR GC metallicities should still be robust. 
If a low-mass galaxy bulge metallicity is underestimated 
due to the presence of younger stellar populations (see 
31]), then the two components of the galaxy may still be 
offset at lower masses. 

Other systems show similar offsets in the aver- 
age MR GC an d bulg e star metallicities. For instance, 
Irlarris fc Harris! (|2002l ) compare the observed metallic- 
ity distributions of MR GCs and resolved bulge stars in 
a massive elliptical galaxy, NGC 5128, with M ste u ar ~ 
6 x 10 10 Mq . The bulge stars in this galaxy at 8 kpc ex- 
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tend to higher metallicities than the MR GCs (see also 
fig. 8 in lBeaslev et alJl2008l ). consistent with the results 
found here. The enr ichment level of Milky Way bulge (e.g. 
IZoccali et al,l I2008T ) is also higher than a typical Galac- 



tic MR GC. This is shown in Fig. [2 
mean metallicity from the ca talogue of 



using a MR GC 



Harris! (1 1 9961) a nd 



Galactic bulge M ste iiar from ICardone fc Serenol ( 20051 ) 



3.2 Discussion 

A difference in the metallicity of two stellar populations 
can reflect a separation in their respective formation 
epochs if the gaseous metals used to form stars were well- 
mixed and the gas enriched over time. Assuming the ob- 
served metallicity difference between galaxy bulge stars 
and MR GCs is a rough timescale indicator, the data 
in Fig. [2] imply that the dominant bulge star-formation 
epoch in massive galaxies was separated in time from the 
peak of the MR GC formation, because the bulges are 
generally more enriched than the typical MR GC they 
host. Since GCs require unique conditions to form, differ- 
ences in the formation epochs of the MR GCs and bulge 
stars likely reflects a change in the star-formation mode 
of the host galaxy. 

The massive galaxy data is consistent with a sce- 
nario where the galaxy experienced a brief epoch of 
high MR GC formation efficiency that was followed by 
a phase where the bulge enriched by ~ 6 0% and grew by 
a sign ficant amount in total M s tellar ■ lHarris fc Harris! 
(2002) speculate a situation like this might occur if GCs 
formation tends to destroy nearby cold-gas clouds (e.g. 
with young high-m ass X-ray binaries found within GCs 
iPower et all l2009h and hence severely limit subsequent 
GC formation. Stellar bulge growth can continue from the 
gas expelled during the peak of GC formation, but GC 
formation will have ended or been severely suppressed 
due to the disruption of large gas clouds. Another pos- 
sibility is that the mode of star formation transitioned 
from a merger-driven star formation mode to one that 
was charac terised by less extreme pressures and/or gas 
turbu lence l)Elmegreen fc Efremov|[l997l : I Ashman fc Zepj 
l200ll ). 

While some intermediate and low-mass galaxies show 
metallicity offsets like those observed in massive galaxies, 
the typical galaxy shows similar bulge star and mean MR 
GC metallicities. If this similarity reflects a coeval forma- 
tion, then bulge star and MR GC formation were roughly 
synchronised. 

In the Milky Way, a relatively large offset in metal- 
licity between the bulge and its MR GC subpopulation 
is found. This might imply the bulge continued to grow 
after the GC formation epoch. 



4 THE RELATIVE FORMATION 

EFFICIENCY OF MR GCS AND BULGE 
STARS 

In this section, the number of MR GCs normalized 
by its host galaxy's Msteiiar observed today is consid- 
ered. This quantity gives some indication of the over- 
all GC formation efficiency in a galaxy, modulo those 
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Figure 3. Relative MR GC numbers (Tmr = 
Ngc /(M a teiiar /10 9 ) against host stellar mass. Ringed 
solid circles have bulge stars that are more enriched by a fac- 
tor of 2 compared to their MR GCs (see Fig. [2j . A substantial 
fraction of the lowest mass galaxies (M ate u ar < 5 X 10 9 Mq) 
contain no detectable MR GC subpopulation. These are not 
shown in this figure, but are given in Fig. [5] 



that are destroyed through dynamical processes (e.g. 
iBaumgardt fc Makindbooj ). 

The relative number of MR GCs in a galaxy, when 
normalized by the host galaxy's stellar mass, is higher 
in massive galaxie s compared to low and intermediate- 
mass galaxies (e. g. iRhode et all2005l ; ISpitler et al.ll2008l ; 
IPeng et al1l2008l ). This quantity is trad itionally referred 
to as the MR GC T-parameter (Tmr', IZepf fc Ashman! 
Il993l ). Fig.[3]shows T M r data for the ACSVCS sample 



The G C numbers were compiled as in ISpitler fc Forbes 
(2009) using the ACSVCS catalogue of IPeng et aL 

(pood) . 

Tmr parameter provides a way to gauge the effi- 
ciency of MR GC production relative to its host galaxy's 
total bulge Mstellar ■ More specifically, it relates GC num- 
bers to the bulge mass accumulated over the entire his- 
tory of the galaxy. To better understand the conditions 
in a galaxy at the peak of GC formation, the relative GC 
numbers normalized by the field stellar mass produced 
during exactly the same epoch should instead be consid- 
ered. This "instantaneous" GC formation efficiency re- 
quires knowledge of the field stars produced at the main 
GC formation epoch, a quantity currently unavailable 
for the ACSVCS galaxies. Thus, only a qualitative as- 
sessment of the instantaneous GC formation efficiency at 
their peak of production can be made. 

For high mass galaxies, the results of Section [3] sug- 
gest that the majority of a galaxy's stellar bulge formed 
in more enriched gas compared to gas used for MR GC 
formation. If GC formation efficiency did not depend on 



© 2010 RAS, MNRAS OOO.riHlTi 



6 Lee R. Spitler 



metallicitjj]] and metallicity is a timescale indicator, then 
the Tmr values for such galaxies are underestimating 
the instantaneous GC formation efficiency. This is be- 
cause the M a teiiar they are normalized against is much 
larger than the M ste u ar produced only with the MR GCs. 
Hence, the instantaneous GC formation efficiency in these 
galaxies will be higher than what their Tmr parameter 
currently portrays. 

These galaxies are highlighted in Fig. [3] The 
M s teiiar in galaxies with no significant temporal bulge 
and MR GC metallicity offset may have formed in con- 
junction with the MR GC subpopulation. Thus, their cur- 
rent Tmr values may provide a reasonable qualitative 
picture of their inherent instantaneous GC formation ef- 
ficiencies relative to massive galaxies. Thus the overall 
effect of moving from the Tmr parameter to an instan- 
taneous GC formation efficiency will largely enhance the 
existing trend: the efficiency of MR GC formation in- 
creases with galaxy Msteiiar ■ 

At the very low-mass end of the galaxy mass distri- 
bution, f5 out of the 25 old ACSVCS galax ies do not hav e 
significant MR GC subpopulations (see iForbea 120051 ). 
although populations of f — 3 MR GCs likely cannot be 
ruled out. This is consistent with the Local Group's 5 
dwarf ellipticals (with Msteiiar ~ 10 8 Mr) that together 
host only two MR GCs l|Sharina et alj|2006h . If the Virgo 
low-mass ellipticals do not show significant MR GC sub- 
populations, then GC production may have been trun- 
cated due to super nova blow-out or r am-pressure strip- 
ping of its gas (e.g. iBoselli et alj|2008h . Low-mass galax- 
ies also generally show high mass-to-light ratios suggest- 
ing they underwent a relatively inefficien t bulge growth 
(e.g iGeha et al.ll2006l ; lBrooks et al.ll2007l 1. thus they pos- 
sibly did not reach the gas densities or star-formation 
rates needed efficiently to produce both bulge stars and 
GCs. 

Alternatively, if f — 3 GCs are actually present in 
some of the low-mass Virgo ellipticals (this cannot be 
ruled out at present), the Tmr values for these galax- 
ies could reach levels sim ilar to those found in massive 
galaxies (|Peng et al ]|2008T ). Whether this implies GC pro- 
duction was enhanced or the normalizing quantity in the 
Tmr parameter, M ste iiar , was underproduced cannot be 
determined without understanding their MR GC subpop- 
ulations, if they even exist. 



5 COSMIC METALLICITY ENRICHMENT 

In this section, the typical formation epoch of the MR 
GCs is estimated using observational constraints on the 
early metallicity enrichment history of galaxies. This will 
enable a comparison to high-redshift galaxy observa- 
tions and further constrain the star-formation histories 
of early-type galaxies. 

A new method is developed to understand the typical 



formation epoch of MR GCs. By comparing the observed 
mean metallicity of MR GCs to empirical age- metallicity 
relationships (AMRs), an estimate of the typical forma- 
tion epoch can be had. GCs provide ideal objects for this 
method of age-dating because they are characterised as 
single stellar populations. In contrast, galaxy observa- 
tions are a complex mixture of multiple stellar popula- 
tions, thus employing an AMR in a similar manner may 
not be so straight-forward. 

Although this method has inherent limitations 
(mostly stemming from the poorly characterised AMRs 
available for redshifts z > 2), it nevertheless provides a 
working approximation of the MR GC formation epoch, 
which can be tested with other observations or simu- 
lations. Also, improved and/or theoretical AMRs (e.g . 



Hemquist & Springcl 



2003; iDave fe Qppenheimerl 12007; 



Kobavashi et al.l 120071 ) can be easily incorporated into 



1 See lForte et al.l (120091 ') for the consequences of assuming GC 
formation efficiency is metallicity-dependcnt. Under this as- 
sumption, they are able to successfully recover a galaxy's ra- 
dial and global stellar populations. 



this analysis, which is relatively simple in design. 



5.1 Analysis 

Three empirical mass-metallicity relationships for star- 
forming galaxy at z = 0.07, 0.7 and 2.2 w ere compiled 
from the lite rature bvlMaiolino et al.l (2008, using the ob- 
servations of Kewlev fe Ellisonll2008l ; ISavaglio et a l. 2005; 
lErb et al ]|2006l ). They also derived emission-line metallic- 
ity measurements in galaxies at z ~ 3.5 and hence effec- 
tively present AMRs for a range of gala xy masses over a 
redshi ft range of z — 0—3.5 or < f 2 Gyrs. Ma nnucci et alj 
(|20QgT ) update the measurement a t z ~ 3.5 with a larger 
sample and use a IChabrierl <|2003l ) initial mass function, 
both of which are adopted in the following. 

Emission-line metallicities from star-forming galaxy 
reflect the dominant enrichment level of the gas being 
used in star formation. Thus by comparing the typical 
MR GC metallicity to these AMRs, the peak redshift of 
their formation can be estimated. Uncertainties, both in- 
trinsic (e.g. possible environmental dependence on the 
AMRs) and systematic (e.g. conversion from observed 
gas [O/H] metallicities to st ellar [m/H] metallicities; see 
iKobulnickv fe Kewlev! |2004| for the conversion adopted 
here, assuming Zq — 0.0f9) translate into non-negligible 
uncertainties on the derived MR GC formation epoch. 
Unfortunately, the extent of these uncertainties are not 
understood, thus the star formation history outlined be- 
low should be considered preliminary. 

Fig. 2] shows the mass-metallicity relationships for 
star forming galaxies at different redshifts. To compare 
the ACSVCS MR GCs to these relationships, the total 
host galaxy M st eiiar o-t the time of their formation must 
be known, since massive galaxies may have continued to 
grow in mass after the peak of MR GC formation (see ij3} . 
As discussed in Section 2] the amount they grew after the 
peak of MR GC formation is currently unconstrained, so 
only a rudimentary correction is attempted to illustrate 
the corresponding implications on the general trends. The 
Msteiiar values of galaxies in Fig. [4] with significant bulge 
star and MR GC metallicity differences (offset by a factor 
of 2, see Fig. [5} have been decreased by a factor of 2. 
One-sided horizontal error bars spanning this factor are 
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Figure 4. Age-dating the formation of MR GC sub-popul ations. Observed em ission-line galaxy mass-metallicity relationships 
(green dashed lines) for the specified redshifts (compiled by Maiolino ct al. 2008). Mean MR GC metallicities are given for Virgo 
early-type galaxies as open circles with a fitted curve shown as a black solid line. Vertical error bars represent the la intrinsic 
metallicity dispersion for the individual MR GC subpopulations, not measurement errors. The solid, large circle represents the 
Milky Way MR GC subpopulation. The relationship between the galaxy-mass metallicity relations (green dashed lines) and the 
fit to the corrected ACSVCS MR GCs (black solid line) suggest the MR GCs formed earlier in more massive galaxies. The Virgo 
galaxy mass-metallicity relationship is given as a red, long-dashed line. Dotted lines are fits to the ±l<r dispersions of the MR 
GC metallicity distributions. Horizontal, one-sided errors show how much certain galaxies moved in the plot, after correcting to a 
^-stellar value that better reflects the M 3te u ar when the MR GCs were forming (see text). 



provided in Fig. [4] for such galaxies and the quadratic fits 
in the Figure reflect this correction. 

In Fig. [4] the fit to the mean MR GC values of mas- 
sive galaxies and the mass-metallicity relationships inter- 
sect at z ~ 3.5, suggesting the peak of MR GC forma- 
tion occurred in massive galaxies around this redshift. In 
lower mass galaxies the intersection implies a peak of for- 
mation between z ~ 2.2 — 3.5. Thus MR GCs typically 
formed in massive galaxies ~ 11.7 Gyrs ago and between 
~ 10.6 — 11.7 Gyrs ago in lower mass galaxies, assuming 
a ACMD cosmology and H = 70 km s" 1 Mpc 1 . The 
lowest mass galaxies do not show a significant MR GC 



subpopulation, hence no constraint is currently available 
from the analysis presented in Fig. [3] 

The intrinsic metallicity dispersions for each MR GC 
subpopulation are also given as vertical error bars in the 
same figure. From the intersection of the mass-metallicity 
relations and the upper dotted line in Fig. [4] the "end" 
of the MR GC formation epoch is z ~ 2.5 for massive 
galaxies and z ~ 1.5 — 2.0 for lower mass galaxies. The 
preferred start of the MR GC formation epoch is z ~ 4—5 
and 2 ~ 3.5 for high and low mass galaxies, respectively. 

As shown in Fig. g] the Milky Way's MR GC sub- 
population has an implied formation age of ~ 12 Gyrs 
(z > 3.5), which is older than similar mass early-type 
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Figure 5. Schematic diagram illustrating bulge formation his- 
tories derived from the analysis of MR GC subpopulations in 
early-type galaxies of the Virgo Cluster. As a rough function 
of the host mass (see Fig. [6jl , these galaxies show three distinct 
bulge formation histories. The "star-formation intensity" is a 
intentionally ambiguous term that is meant to demonstrate 
the fact that the exact mechanism (e.g. local star-formation 
rate, gas turbulence, etc.) that determines the GC formation 
efficiency is unknown. The horizontal dotted line represents 
the threshold for GC production. 



galaxies in the Virgo Cluster. Since the Milky Way's 
bulge appears to be more enriched than its MR GC mean 
the present-day M ste iiar of the bulge may not ex- 
actly reflect its value when MR GCs were forming. In this 
case, the average formation epoch may be closer to z ~ 3. 
This is still older than the cluster early-type galaxies of 
a similar M ate iiar and is contrary to naive expectations 
for chemical enrichment histories, where denser regions of 
the Universe start forming stars at higher redshift. The 
apparent contradiction could reflect the systematic un- 
certainties of employing such preliminary high-redshift 
galaxy mass-metallicity relationships. 



6 DISCUSSION: THE CONNECTION TO 
HIGH-REDSHIFT GALAXY 
OBSERVATIONS 

This section brings together the 3 topics of the preceding 
sections into a preliminary model of bulge star formation 
in Virgo early-type galaxies. Interpretations from direct 
high-redshift galaxy observations are compared to this 
model, followed by caveats and predictions. 

Fig. [5] shows the implied star-formation history of 
early-type galaxies in Virgo as developed from the preced- 
ing analysis on their MR GC subpopulations. The quan- 
tity on the y-axis represents the factor(s) that contribute 
to GC formation. It is intentionally left ambiguous be- 
cause the detailed physics of GC formation are not well- 



understood. Instead the generic term "star-formation in- 
tensity" (SFI) is used as a proxy the mechanisms (e.g. 
star-formation rates, gas turbulence) that dictate GC for- 
mation efficiency within a galaxy. Given on each panel of 
Fig. f5] is a SFI threshold above which GC formation oc- 
curs. 

The scenario in the top panel of the Fig. f5] applies 
to very low-mass galaxies with no significant MR GC 
subpopulation. Here the GC formation epochs are un- 
constrained and the SFI never breaches the threshold for 
GC formation. Detecting strong star bursts in galaxies of 
Msteiiar < 10 9 Af sun at redshifts z < 4 may prove chal- 
lenging. 

In the middle panel of Fig. the SFI history of 
intermediate-mass galaxies is given. These galaxies show 
similar bulge and MR GC metallicities, hence the forma- 
tion epochs of the bulge stars and MR GCs are coeval. 
The peak of SFI in theses galaxies is represented as a 
Gaussian with mean redshift of z ~ 2.5 and spanning 
.z ~ 1.5 — 4 or ~ 1 Gyrs. This age-range matches infer- 
ences made from the observed MR GC metallicity distri- 
bution spread and the observed age-metallicity relation- 
ships (see Fig. 2]). The implied peak of GC and bulge star 
formation apparently coincides with direct observations 
of the cosmic star-formation density peak (at redshifts 
z ~ 2 — 3) , which some believe is dominated by star 
formation in low-mass galaxies (e.g. [Hopkins fc Beacoml 
120061 ; iBouwens et~aT]|2007l ; iReddv et alj|2008h . Here the 
SFI passes the GC formation threshold, to indicate GCs 
formed in these galaxies. 

For massive galaxies with offset bulge star and 
MR GC mean metallicities, the inferred SFI scenario is 
presented in the bottom panel of Fig. [5] The metallicity 
offset implies that the bulge continued to grow after the 
peak of MR GC production, thus a second distribution 
in the SFI history is given at lower redshifts. The 
extent of this bulge formation epoch is unconstrained 
from the present analysis, but is chosen to not overlap 
significantly with the MR GC formation and drop off 
relatively fast to match observations that a significant 



passively evolving 


at z ~ 2 — 3 (e.g. Franx et al. 


2003 


; Ivan Dokkum et al.1 120031; iGlazebrook et al. 


2004 


; Daddi et al. 


2005; van Dokkum et all 20081; 


Cimatti et al.l 2008 


; iBezanson et al.l |2009; Naab et al. 


2009|). Although the SFI was low and few MR GCs 



formed, the bulge likely doubled in total Msteiiar and 
its global metallicity increased by a factor of 2 during 
this second mode of bulge growth. The main GC SFI 
distribution reaches very high values of SFI to signify the 
ultra-efficient period of MR GC production supported by 
the observation (see Fig. f3J). The GC and bulge epochs 
overlap to illustrate the fact that some very enriched 
MR GCs exist. 

In massive galaxies, the peak of MR GC formation 
occurred at z ~ 3 — 4. This redshift is somewhat ear- 
lier than the observed peak in the number density o f 
submillimeter galaxies (z ~ 2.4; IChapman et al.ll2005l ). 
whose properties can be explained b y very intense star- 
formation rates o f ~ 1000 M fl yr" 1 Jllughes et al.lll99Sl ; 



iBlain et al. 1 12002| ; ISwinbank etail |200Sf ). These strong 
star formation events were possibly induced by major 
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galax y mergers llTacconi et al, 20081) or disk insta bilities 



(e.g. lEscala fc Larson! 120081: Ishapiro et all |201(]| ) from 



rapid, cold-gas accretion l|Dekel et al.l 20091 ). In the model 
developed here, this epoch is followed by a mode of 
bulge growth that is not as efficient at GC production. 
Observationally, this period may manifest itself as the 
heavil y dust-obscured massive galaxies foun d at 2 ~ 2 
(e.g. iDev et al.ll2008l : lBussmann et al.ll2009l ). which iust 
finished a subm illimeter-bright phase (see discussion in 
IDev et al.ll200sl ). 

At roughly the same time, approximately z ~ 2.5 or 
~ 11 Gyrs ago, intermediate-mass galaxies were experi- 
encing their peak of bulge growth and MR GC formation. 
The implied SFIs may resemble those in massive galax- 
ies at the same epoch. The average age of bulges in both 
massive and intermediate-mass galaxies should therefore 
be roughly identical, despite a relatively large difference 
in their mean metallicities. 



6.1 Caveats 

The above scenario depends on the following: 

• That there is no significant age difference between 
the galaxy bulge stars and MR GCs. Even though ob- 
viously young galaxies are removed from the analysis, 
transforming galaxy broadband colours into metallicities 
can be complicated by undetected age variations between 
the galaxies and MR GCs. 

• The assumption that the destruction of GCs, be it 
from disrupting tidal forces of the host galaxy or other- 
wise, does not influence the above results significantly. 

• The use of uncertain galaxy age-metallicity relation- 
ships from high-redshift observations. Environmental and 
morphological differences could translate into large sys- 
tematic uncertainties on the attempt to age-date the peak 
of GC formation. 

• The assumption that the apparent offsets in bulge 
star and MR GC metallicities translates into an offset in 
time. Gas accretion, mixing and outflows contribute to 
inhomogeneities in the spatial metallicity distributions of 
a galaxy and will all sabotage the simple model presented 
here. 

The last two points can be better understood by incor- 
porating detailed predictions from theoretic simulations 
of cos mological metal-enrichment and galaxy formation 
(e.g. [Hcrnquist & Sorina'cl 20031: iDave fc Qppenheimerl 
120071 : iKobavashi et al.ll2007i ). These issues will also ben- 
efit from more detaile d observational ana lysis of galaxies 
at high-redshift (e.g. ICalura et ai1l2009h . 



6.2 Predictions 

Presented below are a number of predictions resulting 
from the model outlined above. 

Understanding the mode of bulge star for- 
mation around the transition galaxy stellar bulge 
mass. The three bulge formation histories presented in 
Fig. [5] roughly correlate with galaxy M ste uar , as shown 
in Fig. [6] At M ste iiar ~ 10 10 M the three modes over- 
lap. This bulge M ste u a r corresponds to the characteris- 
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Figure 6. Galaxy M s teiiar histograms for old ACSVCS early- 
type galaxies. Dotted histograms show the total M ate u ar dis- 
tribution of the sample. Solid histograms in upper, middle and 
lower panel correspond to galaxies with no MR GCs, with 
no significant and a significant offset between the typical en- 
richment of its MR GCs and its bulge, respectively. These 
are classified according to the observed metallicity offset be- 
tween their MR GCs and bulge stars, as shown in Fig. [2] At 
Msteiiar ~ 10 10 Mq, the 3 classes overlap. 



tic mass where many g alaxy properties transition (e.g. 
iDekel fc Birnboimll2006l ). This work presents predictions 
for the bulge star formation histories of galaxies below, 
at and above this transition bulge Msteiiar ■ 

GC and galaxy ages. The derived MR GC forma- 
tion redshifts in massive and intermediate-mass galaxies 
suggests an age offset of ~ 1 Gyr should exist in their MR 
GC subpopulations. Typical theoretical modelling and 
measurement uncertainties mean that absolute ages of in- 
dividual extragalactic GCs are difficult to constrain with 
such precision. However, in principle, two very large sam- 
ples of at least 1000 MR GC ages (one for each galaxy- 
mass class) derived from standard Lick absorption-line 
analysis can be used to determine whether or not such 
a relative age offset exists. This minimum sample size is 
derived assuming a mean MR GC age of 11.0 Gyrs (or 
log age = 1.04 in Gyrs) and 11.7 Gyrs (log age = 1.07) in 
intermediate and massive galaxies, respectively. Typical 
uncertainties are assumed to be ±0.15 on t he logarithm of 
the a ge (in Gyrs) of individual GCs (e.g. iProctor et al.l 
120081 ). Note, systematic uncertainties from stellar pop- 
ulation models and the age-metallicity degeneracy may 
complicate such attempts. Furthermore, the spread of 
metallicity may imply a spread of age, thus making this 
prediction more challenging to measure directly. 

Another interesting test of this model is to age-date 
the bulge stars of a massive galaxy (such as Virgo's 
Messier 87) and determine whether the mean age of 1000s 
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of its MR GCs are offset from the bulge's average age by 
~ 1 Gyr. A similar offset is expected in the Milky Way. 

GC formation efficiency. The discovery of an off- 
set between the mean metallicity of MR GCs and their 
host bulge could mean that a significant fraction of the 
bulge M s t e ii ar in massive galaxies formed after their typi- 
cal MR GCs. If the relative MR GC numbers observed in 
galaxies (see Fig. [3| were instead normalized by only the 
galaxy M ste uar that formed in conjunction with the MR 
GCs and it is assumed that GC formation efficiency does 
not depend on metallicity, an extremely MR GC efficient 
formation in massive galaxies is implied. The apparent 
prevalence of ultra-compact dwarfs (which a re sometimes 
consi dered as massive MR GCs; e.g. see iForbes et alJ 
l200ct ) around massive galaxies may be a natural by- 
product of an ultra-efficien t GC formation epoch (see 
discussion in lLarsen|[2009l ). This would mean that the 
local star- formation intensity history, not environment, 
dictates where such mass ive compact star clusters will 
be found l|Hau et alJl200g| ). 
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